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The thermodynamic properties of a spin £ = 1/2 tetrameric Heisenberg antiferromagnetic chain 
with alternating interactions AF1-AF2-AF1-F (AF and F denote the antiferromagnetic and ferro- 
magnetic couplings, respectively) are studied by means of the transfer-matrix renormalization group 
method and Jordan- Wigner transformation. It is found that in the absence of magnetic field, the 
thermodynamic behaviors are closely related to the gapped low-lying excitations, and a novel struc- 
ture with three peaks in the temperature dependence of specific heat is unveiled. In a magnetic 
field, a phase diagram in the temperature-field plane for the couplings satisfying Jaf 1 =Jaf 2 ~Jf is 
obtained, in which various phases are identified. The temperature dependence of thermodynamic 
quantities including the magnetization, susceptibility and specific heat are studied to characterize 
the corresponding phases. It is disclosed that the magnetization has a crossover behavior at low 
temperature in the Luttinger liquid phase, which is shown falling into the same class as that in the 
S=l Haldane chain. In the plateau regime, the thermodynamic behaviors alter at a certain field, 
which results from the crossing of two excitation spectra. By means of the fermion mapping, it is 
uncovered that the system has four spectra from fermion and hole excitations that are responsible 
for the observed thermodynamic behaviors. 

PACS numbers: 75.10.Jm, 75.40.Cx,75.40.Mg 



I. INTRODUCTION 



In recent years, low-dimensional quantum magnets 
have received much attention in condensed matter 
physics. In particular, one-dimensional quantum spin 
chains with competing interactions that show exotic 
physical properties have been extensively studied in the 
past decades. Among others, the dimerized spin-1/2 
antiferromagnetic (AF)-ferromagnetic (F) alternaing 
chain has been widely studied both theoreticallyi* 2 -^ 
and experimentally,— where it was found that the system 
has a gap from the singlet ground state to the triplet 
excited states, and can be mapped onto the 5=1 
Haldane chain 5 ^ if the F couplings dominate. More- 
over, a series of trimerized compounds, including the 
5=1/2 antiferromagnetic F-F-AF chain of 3CuCl 2 -2e!a: 
(ffa;=l,4-dioxane)£ and the ferrimagnetic AF-AF-F 
chains of [Mn(L 2 )(N3) 2 ]n (L=3-methylpyridine) with 
5=5/2, 9 [M(4,4'bipy)(N 3 ) 2 ] n (bipy=bipyridine) with 
M=Co (5=3/2) and Ni (S=l)^ and [Mn(N 3 ) 2 (bpee)] n 
(bpee=trans-l,2-bis(4-pyridyl)ethylene) with 5=5/2, 11 
have been synthesized in experiments. In such trimerized 
spin chains, the topological quantization of magnetiza- 
tion, i.e., the magnetization plateau, has been predicted 
theoretically 12 ' 13 ^ 4 ' 15 ! 16 The predictions were based 
on the theorem proposed by Oshikawa, Yamanaka, 
and Affleckji 7 - which extends the Lieb-Schultz-Mattis 
theorem to give a necessary condition for the appear- 
ance of magnetization plateau in the spin chains with 
translational symmetry. However, the plateau was 
not observed experimently in the above-mentioned 
compounds owing to the weak AF couplings. A plateau 



at m=l/4 (m is the magnetization per site) was later 
observed in an 5=1/2 tetrameric F-F-AF-AF ferrimag- 
net Cu(3-Clpy) 2 (N 3 ) 2 ^^ i 20iSii2 2 . due to the strong AF 
couplings. 

Recently, a spin 5=1/2 tetrameric Heisenberg antifer- 
romagnetic chain (HAFC) with AFi-AF 2 -AFi-F inter- 
actions has been studied, whose ground state was found 
to be in a gapped Haldane-like phase while, importantly, 
it cannot be reduced to an integer spin chaini^ 3 . Twenty 
six years ago, Haldane 24 conjectured that an isotropic 
HAFC with an integer spin has a finite gap from the 
singlet ground state to the triplet excited states, and the 
spin-spin correlation function decays exponentially, while 
the HAFC with half-integer spin has a gapless spectrum 
and a correlation function with a power-law decay. Al- 
though there is no rigorous proof for a general case until 
now, Haldane's scenario has been confirmed experimen- 
tally and numerically in many systems (e.g. Refl25l). Be- 
sides the HAFCs with integer spin, Haldane gap has also 
been found in 5=1/2 spin ladder a 26 i 27 and AF-F alter- 
nating Heisenberg chains* 5 - This is because these spin-1/2 
systems were found to be reducible to an 5=1 HAFC 
when the F couplings are dominantly larger than the AF 
couplings. However, the present 5=1/2 AFi-AF 2 -AFi-F 
tetrameric HAFC cannot be reduced to an integer spin 
chain even if the F coupling dominates, and a gapped 
state with m—0 was observed^ By using a dual trans- 
formation, the Z2 x Z2 hidden symmetry is disclosed to be 
fully broken, and the string order is found non- vanishing 
in the ground state, further suggesting that this spin- 
1/2 tetrameric HAFC system belongs to the Haldane-like 
phase d 8 ' 29 which extends the substance of Haldane's sce- 
nario, namely, the Haldane gap can appear in certain spin 
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half-integer chains. Apart from the gapped state with 
rn=0, a mangetization plateau at m=l/4 was also found 
in this system.— In the critical magnetic fields where the 
magnetization curve is singular, quantum phase tran- 
sitions (QPTs)2£ may happen and consequently, phase 
crossovers are expected at finite temperature. 

As this spin-1/2 AFi-AF 2 -AFi-F tetrameric HAFC ex- 
hibits many interesting behaviors at zero temperature, a 
deeper investigation is still quite necessary. In this pa- 
per, we shall go on elaborating the thermodynamics of 
this system by means of the transfer-matrix renormaliza- 
tion group (TMRG) as well as the Jordan- Wigner (JW) 
transformation, with emphasis on the effects of the cou- 
plings and the external magnetic field on thermodynami- 
cal properties of the system. The possible magnetic phase 
diagram at finite temperature will be presented. The low- 
lying excitations that are closely related to the observed 
thermodynamic behaviors will also be discussed. 

The Hamiltonian of the S—l/2 tetrameric HAFC with 
alternating couplings AF1-AF2-AF1-F in a longitudinal 
magnetic field is given by 

N 

H = y~^(JAFiS4j_3 • S4j_2 + Jaf 2 S4j-2 ' S-tj-i 
j=l 

4AT 

+ JhF 1 ^>ij-i ■ — j^Sij ■ s^j+i) — hyis*,(i) 

3 = 1 

where JaFi, 2 (>0) denote the AF couplings, Jf (>0) de- 
notes the F coupling, and h is the magnetic field. We take 
Jafi as the energy scale and gfis—^ for convenience. iV 
is the total number of the unit cells. The Boltzmann 
constant is taken as fe.e=l. 

The numerical algorithm TMRG method ri^^ which 
is a powerful tool for studying the thermodynamics of 
one-dimensional quantum systems, will be primarily em- 
ployed in the following investigations. As the TMRG 
technique has been discussed in many reviews, we shall 
not repeat the details here. In the following calculations, 
the width of the imaginary time slice is taken as £=0.1, 
and the error caused by the Trotter-Suzuki decomposi- 
tion is less than 10~ 3 . During the TMRG iterations, 
80 states are retained and the temperature is down to 
T=0.02 Jafi in general. In the Luttinger liquid and gap- 
less phases, the temperature is down lower than O.OIJafi 
when calculating the magnetization and susceptibility. 
The truncation error is less than 10~ 4 in all calculations. 

The other parts of the paper are organized as follows. 
In Sec. II, we shall present the TMRG results of the ther- 
modynamic quantities in the absence of magnetic field. 
In Sec. Ill, a magnetic phase diagram at finite tempera- 
ture will be proposed for the case with Jafi= Jaf 2 =Jf, 
and the thermodynamic properties will be discussed in 
the various phases. In Sec. IV, we shall invoke the mean- 
field results from the JW transformation to explain the 
behaviors observed in Sees. II and III. Finally, a sum- 
mary and discussion will be given. 
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FIG. 1: (Color online) Temperature dependence of the specific 
heat of the isolated tetramer systems with (a) Jaf 2 =0; and 
(b) J F =0. 



II. ZERO-FIELD THERMODYNAMIC 
PROPERTIES 

A. Specific heat 

Let us first look at the temperature dependence of 
the specific heat for the isolated tetramer systems with 
Jaf 2 =0 or Jp—0. For Jaf 2 =0, the specific heat decays 
exponentially as T— >0 and has a sharp peak at low tem- 
perature that shifts slightly to lower temperatures with 
increasing Jp, as shown in Fig. 0Ja). For Jf=0, the 
specific heat has a single peak when Jaf 2 /JaFi<1- If 
Jaf 2 exceeds Jaf x , with increasing Jaf 2 the single peak 
splits into double peaks, one of which moves to lower 
temperatures while the other moves to higher tempera- 
tures, as shown in Fig. [lib). The distinct behaviors of 
the specific heat for the two tetramer systems are owing 
to their different energy spectra. When Jaf 2 =0, for each 
tetramer there are four energy levels that decrease with 
increasing Jp, while the gaps among them vary slightly, 
yielding a small shift of the peak. When Jf=0, with 
increasing Jaf 2 the gap between the ground state and 
the first excited states diminishes rapidly, accounting for 
the shift of the low-temperature peak, while the shift of 
the high-temperature peak is owing to the enlarged gaps 
between the energy levels. These limiting cases offer use- 
ful information for better understanding the features of 
the specific heat of the systems with arbitrary couplings, 
which will be discussed below. 

The effects of Jf on the specific heat are firstly dis- 
cussed. In Fig. HJa), the specific heat of the system 
with Jaf 2 / JaFi=1-0 are presented for different Jf- The 
specific heat has a single peak and decays exponentially 
as T— >0, indicating a gapped excitation. With increas- 
ing Jp, the single peak declines and a shoulder appears 
to be gradually prominent at low temperatures. Thus, 
it is expected that the system has at least two gapped 
excitations. With increasing Jp, the gap that is related 
to the lower excitation decreases slowly^ yielding the 
emergence of the shoulder. Figure [Hh) shows the spe- 
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cific heat of the system with Jaf 2 /JaFi=5.0 for different 
Jp. As Jaf 2 i s large, the specific heat has two peaks for 
Jf/JaFi=0.1, like the tetramer system with Jp—0. With 
increasing Jp , a novel peak emerges between the low and 
high temperature peaks, which shifts to higher temper- 
atures with further increase of the F coupling until it is 
merged into the high temperature peak. It can be seen 
that the specific heat behaves rather differently with Jp 
for different Jaf 2 , an( i m certain couplings region, the 
specific heat can exhibit a novel three-peak structure. 

Next, we study the effects of Jaf 2 on the specific 
heat. In Fig. [He), the specific heat of the system 
with Jf/JaFi=1-0 are plotted for various Jaf 2 - When 
Jaf 2 =JaFi, the specific heat has a single peak and de- 
cays exponentially as T^O. When Jaf 2 exceeds Jafi, 
the single peak splits into double peaks, one of which 
shifts to lower temperature side while another moves 
to the higher temperature side with increasing Jaf 2 - 
These behaviors could be understood by means of the 
corresponding tetramer system where the gaps between 
the energy levels are decreasing rapidly with increas- 
ing Jaf 2 <^ It is noticed that the novel peak emerges 
when Jaf 2 is large enough, e.g. Jaf 2 /JaFi=4.0 for 
Jf/JaFi = 1-0. When Jaf 2 continues to increase, the 
novel peak is nearly invariant. When Jp is too large, 
as shown in Fig. [2jd) for Jf/Jafi=10.0, the novel peak 
is absent. These results indicate that the novel peak is 
yielded by increasing Jp after the double peaks have been 
induced by large Jaf 2 • 

The inset of Fig. EJd) shows the parameter regions 
where the specific heat has different peak structures. The 
dashed line separates the regions of the single peak and 
double peaks, while the solid line separates the regions 
of the double and three peaks. The observations in Figs. 
[UJa)-(c) are manifested clearly in this depiction. It is 
shown that when Jaf 2 >JaFi, the single peak starts to 
split into double peaks. When Jaf 2 /Jaf! exceeds about 
3.5, the novel peak can appear with increasing Jf, and it 
would merge into the high temperature peak with further 
increase of the coupling. As the intermediate region is 
enlarged by increasing Jaf 2 , the parameter region of Jp 
for the emergence of the novel peak is wider for larger 
Jaf 2 • It is interesting to point out that the emergence of 
three peaks of the specific heat in the absence of magnetic 
field is nontrivial, which is not the usual feature in low- 
dimensional quantum magnets. 



B. Susceptibility 

The temperature dependence of the susceptibility is 
presented in Fig. [3] for different cases. It is shown that 
the susceptibility has a peak and decreases exponentially 
as T^O. With increasing Jp or Jaf 2 , the peak shifts to 
lower temperatures with the height enhanced, which are 
consistent with the diminution of the gap and the behav- 
iors of the specific heat. As the gap decreases slowly with 
Jf but rapidly with Jaf 2 ^ the susceptibility changes 
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FIG. 2: (Color online) Temperature dependence of the spe- 
cific heat of the tetrameric chain for (a) Jaf 2 /Jaf!=1.0; (b) 
Jaf 2 /Jaf!=5.0; (c) Jf/Jaf^I.0; and (d) Jf/Jaf 1 =W.O 
with JAFa/JAFi^O.5,0.7,0.9, 1.0, 2.0, 5.0, 8.0 and 12.0 from 
top to bottom. The inset of (d) shows the parameter regions 
where the specific heat has different peak structures. 
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FIG. 3: (Color online) Temperature dependence of the spin 
susceptibility of the tetrameric chain for (a) Jaf 2 I Jafi =1-0 
with different Jf; and (b) Jf/Jafi=1-0 with different Jaf 2 - 



slightly with Jp but dramatically with Jaf 2 , as demon- 
strated in Figs. [3]Ja) and (b), respectively. At high tem- 
peratures, the gap is suppressed by thermal fluctuations, 
and the susceptibility goes to coincidence for different 
couplings. These behaviors imply the distinctions of the 
low-lying excitations, which will be reexamined in terms 
of the spinless fermions in Sec. IV. 



III. THERMODYNAMICS IN MAGNETIC 
FIELDS 

In this section, the thermodynamic properties of the 
system with JaFi=Jaf 2 =Jf m the presence of a mag- 
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netic field are studied by means of the TMRG method. 
A phase diagram in the temperature-field plane is pro- 
posed, and the magnetization, susceptibility, and specific 
heat are investigated accordingly in the various phases. 



A. Phase diagram 



As shown in Fig. Ufa), the zero-temperature magneti- 
zation curve is singular at the critical fields h Cl , h C2 , h Ca 
and h s r& suggesting that QPTs^S may happen. These 
transitions arc measured by the divergent peaks of 
dm/dh, which separate the Haldane-like phase, Luttinger 
liquid (LL) phase, plateau phase, gapless phase, and po- 
larized state of the system in a magnetic field. At finite 
temperature, the magnetization plateaus are smeared 
out, and the peaks of dm/dh become analytic, which, 
however, can still describe the crossover behaviors of the 
various phases. Therefore, the magnetization process at 
different temperatures will be studied to obtain the phase 
diagram in the temperature-field plane. 

Figure [Ha) shows that with increasing temperature, 
the two peaks of dm/dh at h Cl (h Cs ) and h C2 (h s ) grad- 
ually merge into a single peak at rather low temper- 
ature, indicating the crossovers from the LL (gapless) 
regime to other regimes. With further increasing tem- 
perature, the high-field peak that separates the plateau 
and spin polarized regimes disappears at T«0.44Jafi, 
and the low-field peak that separates the gapped spin 
liquid and plateau regimes disappears at T«0.68JaFi 
[the inset of Fig. BJa)]. The shifts of the peaks in the 
temperature-field plane compose of the crossover lines. 
Based on the observations, we propose a phase diagram 
in the temperature-field plane as shown in Fig. [4fb), 
from which one may observe that the system has the 
gapped spin liquid, LL, magnetization plateau, gapless, 
spin polarized, and classical phases. 



B. Magnetization 

The temperature dependence of the magnetization 
m(T) in the various regimes are investigated. When 
h<h Cl , the low- lying excitation is gapped. With decreas- 
ing temperature, the magnetization m first increases with 
a power law, then goes down, and finally decays exponen- 
tially to zero with the effective gap A c g=h Cl — h as T^Q. 
As expected, with decreasing the effective gap, the peak 
of m(T) moves to lower temperature side with the am- 
plitude enhanced, as shown in Fig. [5Ja). 

In the LL regime when h Cl <h<h C2 , the gap is closed 
by the field, and the system undergoes a commensurate- 
incommensurate transition at h=h Cl in the ground 
stated At finite temperatures, m(T) shows a minimum 
or maximum at low temperatures, as shown in Fig. [5fb). 
The minima are close to the crossover boundary between 
the gapped spin liquid and LL regimes, while the maxima 
are close to that between the magnetization plateau and 
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FIG. 4: (Color online) (a) Magnetization curves at differ- 
ent temperatures. The inset shows dm/dh as a function 
of h at different temperatures, (b) Phase diagram in the 
temperature-field plane. The solid lines and the dashed line 
are obtained by observing the peaks of dm/dh and the points 
with dx/dT=0, respectively. 



LL regimes, indicating the nonsingular crossovers from 
the LL to the high-temperature regimes. Maeda, Hotta, 
and Oshikawa 38 pointed out that this crossover is univer- 
sal in general gapped one-dimensional (ID) spin systems 
with axial symmetry, resulting from the minus deriva- 
tive of the density of states near the critical field and the 
variation of the Fermi velocity Vp . In Sec. IV, this ob- 
servation would be reproduced by the JW transformation 
and interpreted more physically. 

In Ref. [H, a linear dependence of the transition tem- 
perature (T m ) on the field T m =xo(h — h c ) (h c is the criti- 
cal field where the gap is closed) near the critical field h c 
was proposed for the S=l Haldane chain. The field de- 
pendence of the crossover temperature for the present 
tetrameric HAFC is shown in the lower inset of Fig. 
03L). It can be seen that the crossover temperature T m 
varies like a sine function of the field and behaves lin- 
early near the critical fields /i C] and h C2 . The coefficient 
in— 0.76238 proposed in Ref. |38| for the S=l Haldane 
chain near h c fits well to our data near both h Cl and 
h C2 (see the solid lines in the inset), indicating that this 
tetrameric chain slightly above h Cl and below h C2 could 
also be well described by the free fermion theor y 39 i 40 i 41 
and the crossover in this LL phase is of the same class 
as that in the Haldane chain. Although this tetrameric 
chain cannot be reduced to a typical Haldane chain with 
an integer spin, their analogous dispersion relation of the 
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low-lying excitations are revealed by the same class of 
this crossover behavior of magnetization. 

In the magnetization plateau, m approaches 0.25 as 
T^O. As shown in Fig. El[c), when h C2 <h<h m ~ 
I.IJaFi, fn increases slowly with cooling temperature to 
a certain value and, then, rises rapidly to 0.25, where 
h m is the crossing field of the magnetization curves at 
low temperatures in the plateau states, as marked by the 
triangle in Fig. 01a). When h m <h<h C3 , m increases 
to a maximum and then declines to 0.25. The different 
behaviors are also observed in the magnetization curves 
at various temperatures [Fig. Ufa)]. In the field h m , m 
changes slowly at low temperature, like the m(T) curve 
with /i/Jafi=1-1 in Fig. 03c). These magnetic behaviors 
in the plateau state have also been noted in the spin-1/2 
trimerizecU^ and F-F-AF-AF tetrameric chains^ In Sec. 
IV, it would be found that these common features in the 
plateau states may result from the crossing of fermion 
and hole excitations. The field h rn is in the middle of 
the plateau and corresponds to the midpoint of the gap 
between the crossing fermion and hole spectra. When 
h>h m , the maximum of m(T) with dm/dT=0 measures 
the change of the magnetic properties. For h<h m , al- 
though m keeps declining, such a change also exists, 
which can be visible in the temperature dependence of 
dm/dT. Corresponding to dm/dT=0 for h>h m , dm/dT 
has maxima that are marked by solid squares to separate 
the high temperature regime for h<h m , as shown in the 
inset of Fig. [5jc). As the field approaches h m , the max- 
imum tends to disappear. 

In the gapless regime, the gap in the plateau state can 
be closed by increasing the field. Thus, the magnetic be- 
havior m{T) is analogous to that in the LL phase, where 
there are also minimum and maximum observed at rather 
low temperatures [Fig. [HJd)]. Such a crossover in the 
gapless phase between a m^O plateau and the saturated 
state has not been reported. The magnetic behaviors 
of the spin-1/2 F-F-AF antiferromagnetic chaini^ have 
been studied numerically, but no such a crossover is ob- 
served between the m=l/6 plateau and the saturated 
state when the temperature is down to 0.025Jf. It is 
also noticed that the crossover temperatures in the gap- 
less regime are lower than those in the LL phase, which 
would be interpreted in terms of the spinless fermion in 
Sec. IV. 



C. Susceptibility 

The behaviors of the susceptibility x hi the various 
regimes will be discussed in this subsection. When 
h/ JaFi is less than about 0.3 in the gapped spin liquid, 
X has a single peak at T~A (A is the gap in the absence 
of magnetic field), and approaches zero exponentially as 
T—>0. With further increasing h, a new peak emerges and 
moves to lower temperatures with the height enhanced, 
while the high temperature peak becomes smoother, as 
shown in Fig. [f^a)- In the high temperature region T>A, 
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FIG. 5: (Color online) Magnetization as a function of temper- 
ature for (a) h<h ci ; (b) h cl <h<h C2 \ (c) h C2 <h<h C3 ; and (d) 
h C3 <h<h s . The insets show (b) the minima (maxima) and 
the field dependence of the crossover temperature; (c) the be- 
havior of dm/dT for h/J AFl =0.9, 1.0, 1.02, 1.04, 1.05 and 1.1 
from bottom to top; and (d) the minima and maxima in the 
gapless regime. 



the susceptibility under different fields coincide because 
the gap is suppressed by thermal fluctuations. In the 
fermion mapping in Sec. IV, the system has two positive 
energy excitations (fermion excitations) and two minus 
energy excitations (hole excitations) when h=0. With 
the decrease of the excitation energies induced by in- 
creasing the field, the susceptibility contributed from the 
fermion excitations moves to lower temperatures with the 
amplitude enhanced, while that from the hole excitations 
shifts to higher temperatures with the height decreased, 
both of which are responsible for the behaviors observed 
in Fig. [UJa). In the LL, x is finite as T— >0, as shown in 
Fig. [DJb). With cooling temperature, x increases slowly 
to the temperature T~A and then has a sharp rise until 
to a peak at rather low temperatures, which results from 
the closure of the gap. 

In the plateau state, due to the open of a gap, the 
susceptibility has a single peak, and approaches zero ex- 
ponentially as T^0. With increasing the field, the peak 
moves to higher temperatures with the height declined 
when h<h m , while it shifts to lower temperatures with 
the height enhanced after h exceeds h m , as shown in the 
inset of Fig. [|3a). The field dependence of the peak tem- 
perature is displayed by the dashed line in Fig. [4jb), 
where the field with the highest peak temperature is h s , 
as marked by a triangle. These distinct behaviors of the 
susceptibility also result from the crossing of the fermion 
and hole spectra. When h<h m , x IS dominated by a 
hole branch whose gap is enhanced with increasing the 
field, yielding the peak to move to higher temperatures 
with the height decreased. After h exceeds h m , the two 
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FIG. 6: (Color online) Temperature dependence of the sus- 
ceptibility for (a) /i/Jaf!=0.1, 0.2, 0.3, 0.34, 0.36, 0.38, 0.4 and 
0.42 from bottom to top; (b) /i/Jaf^O.5 and 0.55; the inset 
of (a) h C3 >h>h C2 ; and (b) h s >h>h C3 . 
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spectra cross and the susceptibility is dominated by the 
fermion branch whose gap declines with the increasing 
field, yielding the peak to move to lower temperatures 
with the height enhanced. Different from the gapped spin 
liquid, the susceptibility does not show double peaks in 
the plateau phase. In the gapless regime, owing to the 
closure of the gap, the susceptibility exhibits the same 
features as those in the LL, as shown in the inset of Fig. 

Ep>). 



D. Specific heat 

In this subsection, the specific heat is explored in de- 
tail. When h<h Cl , the specific heat C has a single peak 
and approaches zero exponentially as T^0. The peak 
shifts to higher temperatures with the height decreased 
when A ff diminishes [Fig. Efa)]. This shift is attributed 
to the hole excitations whose gaps increase with increas- 
ing the field. It is noticed that the peak temperature, 
Tp e ak/ Jaf ± — A, in this system is distinct from the re- 
sult Tp ea k/J— 2A' of the 5=1 Haldane chain (J and A' 
are the coupling and gap of the 5=1 Haldane chain, 
respectively) M- As the field h approaches h Cl , a shoul- 
der gradually emerges at low temperature, which is a 
signature of approaching the quantum critical poin t 30 i 36 
and is from the hole excitations. 

In the LL, the linear temperature dependence of the 
specific heat at low temperature is observed. With a fur- 
ther increase of the field, the system shifts away from 
the quantum critical point h Cl . Thus, the shoulder at 
low temperature is smoothed down gradually, as shown 
in Fig. \T[b). The disappearance of the shoulder is at- 
tributed to the hole excitations. 

In the plateau state, the specific heat decays exponen- 
tially when T^0 because of the open of a gap. Near 
the lower critical field h C2 , the shoulder that emerges in 
the LL vanishes gradually, while near the upper criti- 
cal field h C3 , a double-peak structure emerges. As ex- 
pected, the crossing of the fermion and hole spectra af- 



FIG. 7: (Color online) Specific heat of the system with (a) 
h ci >h; (b) hc2>h>h c i; (c) h c:i >h>h C 2; and (d) h s >h>h C 3- 



fects the behavior of the specific heat. With increasing 
the field, the peak of the specific heat moves to lower 
temperatures when h<h mi and when h exceeds h m , the 
peak starts to move to higher temperatures, as shown 
in Fig. EJc). However, different from the magnetiza- 
tion and susceptibility, the behavior of the specific heat 
cannot be characterized simply only by the hole (h<h m ) 
or fermion (h>h m ) excitations, although the crossing in- 
deed changes the features of the behavior. This is because 
the magnetization and susceptibility are determined only 
by the occupied number of the excitations but the spe- 
cific heat is affected by both the numbers and energies of 
the quasiparticles. 

In the gapless regime, with increasing the field, the 
high temperature peak of the specific heat keeps nearly 
intact, while the low-temperature peak that occurs near 
the critical field h C3 in the plateau state moves to lower 
temperatures with the height declined [Fig. |7[d)] . These 
behaviors will be analyzed in the next section. 



IV. JORDAN- WIGNER TRANSFORMATION 
AND SPINLESS FERMION MAPPING 

In order to explain the thermodynamic behaviors ob- 
served in the above sections, the elementary excitations 
of the system are studied using the JW transformation. 42 
The 5=1/2 spin operators can be transformed into the 
spinless fermion operators through JW transformation 



5+ = cj e " r £,<* c K St 



(4c - \), (2) 



where cj and Ci are the creation and annihilation opera- 
tors of the spinless fermion, respectively. For this 5=1/2 
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tetrameric HAFC, four kinds of spinless fermions should 
be introduced: 



bib 



m<j 



'43-1 



■n 



43-3 



a' a 



3 "3 



43-2— °3°3 2 ' 
It, 



C*2 — r'r 1 Z — d' d 

°4j-l ~ S 3 2' °4j — a j a 3 2' 



(3) 



After the JW transformation, the interactions of 
the original Hamiltonian are transformed to the nearest- 
neighbor hoppings of the fermions, and the Ising terms 
become the nearest-neighbor density-density interactions 
that will be treated by the Hartree-Fock (HF) approxima- 
tion. By performing a cumbersome derivation, we obtain 
a mean-field Hamiltonian after omitting the constant 



N 



H H f = X]^[ jAF i( db ~ h ~ M d <i~ \)}a)aj 

3=1 

+ 2[ j a Fi K- -) + JAF*{dc--£]b\bj 

+ \[JAF 2 {db~ -) + JhFMd- 



J t^i(\ -PAB)a]bj + Jaf 2 (^ ~PBc)b]cj 
+ JafAt; -PcD)c\dj - J F (- - p D A)d\a j+1 



2 ty 2 i JUA l u, j u, 3' 

N 

+ h.r.) hY.'"]"- ■ h ?j ■ 'Yj ■ ''/I,:. (4) 

3=1 

where the occupied fermion numbers are d a =(ajCij) , 
db=(bjbj), d c ={CjCj), dd=(djdj), and the covalent bond- 
ings are pAB=(b]aj), pBC=(c]bj), p CD ^(d f J c j ), and 
PDA={a] + \dj). The brackets (• • •) denote either the 
HF ground state average (T=0) or the thermal average 
(Ty^O). By making the Fourier transform and then Bo- 
goliubov transformations, a quadratic Hamiltonian can 
be obtained 



where a, (3, 7, and A denote four excitation spectra with 
the dispersion relations Lo k (i=a, (3, 7, A). 



A. Spin gap and magnetization 

In the mean-field calculations of the ground-state prop- 
erties, the occupation numbers and covalent bondings 
are self-consistently calculated by minimizing the ground 

state energy with the constraint (S , f t) = (X)i=i <Sf) = 0; 
i.e., the total number of spinless fermions is 2N. The self- 
consistent calculations give rise to four excitation spectra 



[Fig. HJa)]. In the absence of magnetic field, ui k - 



and 
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u>2- Thus, the ground state is obtained by fill- 



ing up the two negative spectra a and (3, The gap from 
S* ot =0 to Sf ot =±l, corresponds to the energy of adding 
or removing a fermion 

A ff F = ^ = ^{A-{A 2 -4[Ji Fl (l + 2p AB f 

x (l + 2p 2 CD ) -2J| Pi Jap 2 Jf(1 + 2p ab ) 

x (1 + 2 Pb c)(1 + 2 P cd){1 + 2pda) 

+ J 2 AF2 4(l + 2PBc) 2 (l + 2pDA) 2 ]}hK (6) 



where the parameter A is 



A = J\ 
+ J 2 



XFl [(l + 2p AB ) 2 
(1 + 2p BC ) 2 



AF 



- (1 + 2p CD 

jf(l + 2p 



DA 



(J) 



For Jafi =Jaf 2 =Jf=1, the mean-field result gives the 
gap A//i?=0.4377JAFi, which agrees with the DMRG re- 
sult A=0.435J AFl ^ 

In Fig. [H^b), the gap is plotted as a function of 
Jafi • The HF results agree with the DMRG values when 
JAFi/ Jf>0-7. In Figs. EJc) and (d), Jf and Jaf 2 depen- 
dences of the gap are plotted, respectively. It can be seen 
that when Jp and Jaf 2 are smaller than Jafi > the HF 
results are consistent with the DMRG results, but with 
increasing Jp or Jaf 2 j the HF results become worse when 
Jf or Jaf 2 is prominent. 

In a magnetic field, new quasiparticles can be excited 
in the 7 and A branches. It is found that the excited 
fermions can interpret the magnetization process at zero 
temperature. When the field is less than the gap, no 
fermion is excited and m—0. When the field closes the 
gap at h Cl , new fermions are excited in the 7 branch, 
and m increases until the 7 branch is fully filled at h C2 . 
The gap between the spectra 7 and A is the width of the 
plateau, say, h Ca — h C2 . When the field exceeds h C3 , the 
fermions in the spectrum A are excited, and the system 
is fully spin-polarized. 



B. Zero-field thermodynamics 

At finite temperature, the numbers of the excitations 
obey the Fermi distribution function n k = l/(e Wfc / T + 1) 
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J /J 
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FIG. 8: (Color online) (a) The elementary excitation spectra 
of the system. The gap as a function of (b) Jaf ± ; (c) Jf; and 
(d) J A f 2 . The DMRG results are from Ref. [23|. 



(i = a,/3, 7, A), from which the thermodynamics can be 
obtained by the self-consistent calculations for the occu- 
pation numbers and covalent bondings. 

In Fig. [HI the specific heat obtained from the HF cal- 
culations are displayed, which qualitatively agree with 
the TMRG results. With increasing Jp or Jaf 2 , the gap 
is decreased. Therefore, the peaks of the specific heat 
contributed from and 7 spectra move to lower temper- 
atures with the height decreased, while those from a and 
A spectra shift to higher temperatures with the height 
declined, both of which compose the observed shoulder 
and double peaks in the specific heat. As the spectra f3 
and 7 shift slowly with Jp but rapidly with Jaf 2 i the 
specific heat exhibits a shoulder with increasing Jp but 
double peaks with increasing Jaf 2 a t low temperature. 
The three peaks in the specific heat that cannot be re- 
produced in the mean-field theory may result from the 
interactions of the quasiparticles, which might induce ex- 
citations between a(X) and /3(7). It is expected that the 
energies of the induced excitations are enhanced with in- 
creasing Jf , yielding the observed shift of the novel peak 
from the low to high temperatures in Fig. [DJb). 



C. Thermodynamics in magnetic fields 

At zero temperature, the spectra a and (3 are fully oc- 
cupied, while 7 and A spectra are empty in the absence 
of magnetic field. Thus, the excitations at finite temper- 
ature are of the hole (a and /?) and fermion (7 and A) 
types, respectively. In a magnetic field, the energies of 
the spectra decrease, yielding the excitations with minus 
energies to become hole excitations at finite temperature. 




FIG. 9: (Color online) Temperature dependence of the spe- 
cific heat for (a) Jf/Jafi =0.5, 1.5, 2.5 and 3.0; and (b) 
Jaf 2 /^af 1 =0.5, 1.0, 1.5 and 2.0 from top to bottom. The 
comparisons with the TMRG results are shown in the insets. 



In this subsection, some typical thermodynamic behav- 
iors could be explained by the combination of the hole 
and fermion excitations. 

When h Cl <h<h C2 and h C3 <h<h s , the minimum and 
maximum of m(T) at low temperature are observed in 
the mean-field results, and the crossover temperatures 
in the gapless regime are lower than those in the LL, 
as shown in Fig. fTUT a). As the magnetization is pro- 
portional to the total number of the fermions, the min- 
imum and maximum indicate that the hole and fermion 
excitations dominate at low temperature, respectively. 
In the vicinity of the lower critical field h Cl (h C3 ), the 
fermion spectrum 7 (A) crosses the Fermi level slightly, 
as shown in the lower inset of Fig. fTUT a). making the 
excitations with minus energies become hole type. Thus, 
the low-temperature behavior of m(T) is determined by 
the competitions of the fermion and hole excitations near 
the Fermi level. Due to the dispersion relation, the den- 
sity of states D(uj) of the holes is larger than that of the 
fermions near the Fermi level, yielding more holes excited 
and thus the decrease of m at low temperature. After the 
few holes are occupied, the fermion excitations dominate 
and, thus, a minimum of m emerges. The analogous ar- 
guments for the maximum near the upper critical field h C2 
(h s ) also apply. The lower crossover temperature in the 
gapless regime is attributed to the larger values of the 
gradient of the Fermi velocity with respect to the field 
dvF/dh of the spectrum A near the Fermi level, which 
diminishes the difference of the density of states in the 
vicinity of the Fermi level. 

Different from the magnetization, the fermion and hole 
excitations have equivalent contributions to the suscepti- 
bility. However, with increasing the field, the fermion ex- 
citations move close to the Fermi surface, while the holes 
move away. Thus, the peak of the susceptibility from 
the excited fermions shifts to lower temperatures with 
the height enhanced, while that from the excited holes 
exhibits opposite behaviors, as shown in Fig. [TWb). All 
the observed behaviors of the susceptibility in the gapped 
spin liquid and plateau phases can be understood on the 
basis of the above analyses, as mentioned in Sec. III. 
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FIG. 10: (Color online) Mean-field results in different mag- 
netic fields for (a) magnetization m(T) at low temperatures 
in the gap less regime; (b) susceptibility contributed from the 
fermion and hole excitations. The lower inset of (a) depicts 
the dispersion in the vicinity of the lower critical field. 



Now let us discuss breifly the specific heat. The spe- 
cific heat can be divided into those from the fermion and 
hole excitations, which are denoted as C/ and Ch, re- 
spectively. It is found that C/ and Ch behave rather dif- 
ferently with the field, and Ch contributes essentially to 
the observed characteristic features, which are not shown 
here for brevity. For example, when h<h Cl , the gaps of 
the hole excitations a and /3 are enhanced with increasing 
the field, making the peaks of Ch move to higher tem- 
peratures with the amplitude decreased. As h approaches 
h Cl , a low temperature peak in Ch emerges, yielding the 
shoulder in the specific heat. In other phases, similar 
analyses are also applicable. 

As discussed above, the spinless fermion mapping can 
be used to explain both the ground-state and thermo- 
dynamic properties in a wide range of the parameters. 
The fermion and hole excitations provide a simple way 
to understand the complex thermodynamic behaviors of 
the system. 



V. SUMMARY AND DISCUSSION 

The thermodynamic properties of the spin 5=1/2 
tetrameric HAFC with alternating couplings AF1-AF2- 
AFi-F have been studied by means of the TMRG method 
and JW transformation. In the absence of magnetic 
field, the thermodynamic behaviors are determined by 
the gapped low- lying excitations. The specific heat can 
have single peak, double peaks, and three peaks for differ- 
ent couplings. The shoulder and double peaks in the spe- 
cific heat are attributed to the decreased energies of the 
gapped excitations. The novel intermediate peak in the 
specific heat results from the increase of Jf after the dou- 
ble peaks have been induced by large Jaf 2 • With further 
increase of Jf, the novel peak shifts to higher tempera- 
tures and finally merges into the high-temperature peak. 
The susceptibility is found to have a peak that shifts to 
lower temperatures with the amplitude enhanced with 
increasing Jf or Jaf 2 - 



A phase diagram in the temperature-field plane for 
JaFi=Jaf 2 =Jf is obtained by analyzing the crossover 
behaviors of the various phases at finite temperature. 
The system is unveiled to contain the gapped spin liq- 
uid, LL, magnetization plateau, gapless, spin polarized, 
and classical phases. 

In the LL and gapless regimes, the magnetization curve 
exhibits a minimum or maximum at low temperature, 
representing a nonsingular crossover. The linear field de- 
pendence of the crossover temperature in the LL with the 
same ratio as that in the S=l Haldane chain indicates 
that the crossovers in the two systems belong to the same 
class, which implies that the Haldane-like phase of this 
tetrameric HAFC has an analogous low-lying dispersion 
relation to that in the 5 = 1 Haldane chain. In the plateau 
state, as a crossing of a fermion and a hole spectra at 
h m , the temperature dependence of the magnetization 
behaves differently below and above the field. 

When h<h Cl , the susceptibility has double peaks with 
increasing the field. In the plateau state, the susceptibil- 
ity has a single peak, and owing to the crossing of the 
spectra, the peak moves to higher temperatures when 
h<h m , and to lower temperatures when h>h m with in- 
creasing the field. In the LL and gapless regimes, \ is 
finite as T^0. The susceptibility in any fields coincides 
when T>A due to the thermal fluctuations. 

In the gapped spin liquid, the specific heat has a shoul- 
der as the field approaches the critical field h Cl . In the 
LL, the specific heat as a function of temperature be- 
haves linearly at low temperature, and the shoulder is 
smoother down gradually by increasing the field. In the 
plateau, owing to the crossing of spectra, the peak of the 
specific heat moves to lower temperatures when h<h m , 
and to higher temperatures when h>h m with increas- 
ing the field. As the field approaches h C3 , a small peak 
emerges at low temperatures. In the gapless phase, with 



increased field h 



the low-temperature peak moves 



to lower temperature with the amplitude decreased. 

By means of the JW transformation and mean-field 
approximation, the low-lying excitations and thermody- 
namic properties of the system are studied to understand 
the TMRG results. It is unveiled that the system has four 
excitation spectra with a gap, which may account for 
the magnetization process at zero temperature. At finite 
temperature, it is found that the thermodynamic behav- 
iors are determined by the combination of the fermion 
and hole excitations. The complex thermodynamic be- 
haviors in the TMRG results can be understood well 
within the free fermion mapping. 

Finally, we would like to add that the observations pre- 
sented in this paper for this spin-1/2 tetrameric HAFC 
would be expected to test experimentally in future. 
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